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Abstract 

The effect of symmetry energy and cross section had seen on the fragment production in the mul- 
tifragmentation of ^''A^eio+^'^ A^cio ^-iid ^^"^ Auyg ^^^"^ Aujg at incident energy 50-1000 MeV/nucleon 
using isospin dependent quantum molecular dynamics model. To see the effect of symmetry en- 
ergy and cross-section in more effective way relative yield of ^^'^ Aw/g and ^"^A^eio was studied. It is 
observed that relative yield of free nucleon, light mass fragment at isospin dependent cross-section 
is almost same at different symmetry energy values. While for fixed cross-section relative yield 
is influenced by symmetry energy. We compare the stiffness of symmetry energy suggested by 
different groups with our value. To verify our result we compare our study with ALADIN data. 

PACS numbers: 25.70.-z, 25.70.Pq, 21.65.Ef 



* Electronic address: 



suneel.kumar@thapar.edu 



I. INTRODUCTION 

The science of nuclear physics deals with the properties of nuclear matter. The study of 
nature of matter and strength of nuclear interaction is key to understand many fundamental 
problems. Understanding the behaviour of nuclear matter at density and isospin away from 
normal nuclear matter (T ^ MeV, po ~ 0.16 fm^^, N ^ Z ) has gained tremendous im- 
portance In the past decade only light ions or particles could be accelerated to produce 
the collision. But it is possible to accelerate the heavy ions to a large amount of energy. 
The term heavy ion is generally used for nuclei which are heavier than the helium nucleus. 
The heavy ion collisions has attracted the scientific world because of its various features. 
The intermediate energy heavy ion collisions are helpful in studying the nuclear matter at 
extreme conditions of temperature and pressure, which can be correlated with astrophysical 
happening like supernova explosion, neutron star etc. 

Multifragmentation is one of the important phenomena occurring at intermediate energies. 
The process of breaking of colliding nuclei into several small, medium size fragments is 
called multifragmentation. Those fragments which are highly excited and neutron-rich sub- 
sequently undergo de-excitation to cold and stable isotopes. Similar hot nuclei are also 
produced in the interior of a collapsing star and subsequent supernova explosion [2|. 
The production of these nuclei depends on thejr internal excitation and^is sensitive to the 



symmetry energy part of the binding energy [3]. Puri and co-workers 



^ has analysed the 



time evolution of fragment production for symmetric nuclei Ca + Ca, Xe + Xe and Au + 
Au. There are many parameters which effects the multifragmentation. Symmetry energy is 
one of the most promising parameter to study the multifragmentation after the collision be- 
tween two nuclei. Symmetry energy Egymip) of nuclear matter characterizes how the energy 
changes as one move away from equal numbers of neutrons and protons. In the heavy ion 
collisions the dynamics of the collisions between any two nuclei is also sensitive to isospin 
dependent NN cross section. So, it is an interesting and important goal of heavy ion physics 
to extract information of symmetry energy and its density dependence. Many groups had 
found the relation of symmetry energy [5|-[7|. They consider different values of E^y^ and 7 
and see the influence of stiffness and softness of the symmetry energy. The best estimate of 
the density dependence of the symmetry energy can be parameterized as ref. 8| 

Esym = E )''' (1) 

Po 

In multifragmentation the measurement of fragment isotopic yield distribution can provide 
important insight into the symmetry energy and the decay characteristics of these nuclei. It 
has been shown in experimental measurements that the ratio of the fragment yields, Ru (N, 
Z), taken from two different multifragmentation reactions, 1 and 2, obeys an exponential 
dependence on the neutron number (N) and the proton number (Z) of the fragments; an 
observation known as isoscaling |9i]. The dependence is characterized by the relation 

i?2i(iV, Z) = Y^iN, Z)/Y,{N, Z) = Ce("W) (2) 

Where Y2 and Yi are the fragment yields from the neutron-rich and the neutron-deficient 
systems, respectively, C is an overall normalization factor, and a and (5 are the parameters 



characterizing the isoscahng behaviour. So far, the isoscahng behaviour has been studied 



experimentally and theoretically for different reaction mechanisms 1CH12|. D. V. Shetty et 



ai, 



13 



14| , experimentally understand the correlation between the temperature, density 



and symmetry energy of multifragmentating system as it evolves with excitation energy. 
They also study the relative reduced neutron and proton densities as a function of excited 
energy of fragmenting source for the ^^Fe + ^^Ni and ^^Fe + ^^Fe reaction. Their study 
shows a steady decrease in reduced neutron density and an increase in proton density 
with increase excitation energy. Two nucleons can collide if they come closer than certain 



distance. Multifragmentation is also influenced by cross section. The nuclear cross-section 



is used to characterize the probability that a NN collision will occur. C. Zeitlin et ai, 15 1 



experimentally study the effect of cross section on fragment production. R. Wada and 



co-worker 16| shows that value of cross-section decrease with increase in incident energy. 



In this paper we compare the various forms of symmetry energy suggested by different 
groups with our suggested value and see its softness and stiffness for different values of the 
parameter 7. We have tried to understand the influence of symmetry energy on multiplicity 
of fragments for neutron rich and deficient nuclei at different incident energies i.e. we see 
the effect of mass dependence on fragment production under the effect of symmetry energy. 
In addition to symmetry energy we analyse the relative production of fragment FN (A 
= 1) and LMF (2 < A < 4)) under the effect of cross section ( isospin dependent and 



fixed (20 and 50 mb) cross section) . To verify the result we compare our calculations with 



experimental ALADIN data 



18| . This study is done within the frame work of isospin 



4 



quantum molecular dynamics model 



191]. Section II deals with the model, Sec. Ill discusses 



the results. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 



Our study is performed within the framework of IQMD 19|] model, which is an improved 



version of the QMD model 



20| developed by J. Aichelin and coworkers and then improved 



by Puri and coworkers has been applied to explain various phenomenon such as collective 
flow, disappearance of flow, fragmentation and elliptical flow successfully. The isospin degree 



of freedom enters into the calculations via symmetry potential, cross-sections and Coulomb 



interactions |19| . The details about the elastic and inelastic cross-sections for proton-proton 



and neutron-neutron collisions can be found in ref. 



19| . In IQMD model, the nucleons 



of target and projectile interact via two and three-body Skyrme forces, Yukawa potential 
and Coulomb interactions. In addition to the use of explicit charge states of all baryons 
and mesons, a symmetry potential between protons and neutrons corresponding to the 
Bethe- Weizsacker mass formula has been included. Skyrme forces are very successful in 



the analysis of low energy phenomena such as fusion, fission and cluster-radioactivity, where 



nuclear potential plays an important role 21 1. 



where hadrons propagate with Hamilton equations of motion: 



drj 
'dt 



d{ H ) 
dpi 



dpi 
dt 



d{ H ) 
dvi 



(3) 



5 



with 



H) = {T)+{V) 



X fj {r', 0, t)drdr' dpdf/ 



(4) 



The baryon-baryon potential in the above relation, reads as: 



Skyrme ' Yukawa 



' ^Coul ^ ^ XII- 



sym 



ti5{f' - f) + t25{f' - r)p^'^ 



r' + r 



+ t; 



exp{\r' — r|//i) 
(|P-r1//i) 



I ly^ 



+t6-T3%^5(rt 



^'0 



(5) 



Here and Zj denote the charges of i*^ and j*'^ baryon, and TJ, T| are their respective 



th 



T3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 
Coulomb interaction only. The binary nucleon-nucleon collisions are included by employing 
the collision term of well known VUU-BUU equation. During the propagation, two nucleons 
are supposed to suffer a binary collision if the distance between their centroids 



\ri -'rj\< xl^jCTtot = cr{y/s,type) 



(6) 



"type" denotes the ingoing collision partners (N-N, N-A, N-tt,..). In addition, Pauli blocking 
(of the final state) of baryons is taken into account by checking the phase space densities 



in the final states. The phase space generated using IQMD model has been analysed using 



minimum spanning tree algorithm 21 1 
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FIG. 1: Schematic diagram to show the stiffness of symmetry energy suggested by various groups. 
III. RESULTS AND DISCUSSION 



We have simulated 1000 events involving the symmetric reactions of '^^Nciq A^Cio and 
^^Mmtq +19^ Au79 at incident energies 50, 100, 200, 400, 600, 1000 MeV/nucleon for central 
collision (scaled impact parameter is 0.3) using soft equation of state. Here, we use three 
different forms of symmetry energies. For all the calculations we take symmetry energy 
corresponding to normal density is 32 MeV and value of 7 which characterizes the stiffness 
of the symmetry energy is 0, 0.66 and 2. In figure 1, we compare the stiffness of symmetry 
energy suggested by different groups. 

Red and royal coloured symbol having 7 = 2 are highly stiffed than others. Here we anal- 
yse that larger the 7 value more the stiffness. For figure 1 density dependence of symmetry 
energy obtained from various groups shows the close agreement with parameterized forms 

7 
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FIG. 2: Influence of symmetry energy on fragment production for symmetric reaction Au + Au 
and Ne + Ne at different energies for central collision ( 6=0.3). 



of the density dependence of the symmetry energ y, g iven as in eq. (1), where ^ 31-32 



MeV and 7 ~ 0.55-0.69 which is agree with ref. 



13| . We study the influence of symmetry 



energy on fragment production at different incident energies in figure 2. On the basis of 
mass number the fragments are classified as free nucleons [ A = 1 ], light mass fragments 
(LMFs) [2 < A < 4], intermediate mass fragments (IMFs) [5 < A < Atot / (3 or 6)]. 
From the Figure 2, we observed that the multiplicity of FN ) LMF ) IMF for both neutron 
rich and deficient nuclei. On considering the mass effect for lightest nuclei (Ne -|- Ne) at low 
energy (50 MeV/nucleon) the relative ratio of FN: LMF: IMF is 8: 3: 1 and for high energy 
(1000 MeV/nucleon) is FN: LMF: IMF :: 7: 2: 1. Due to less Coulomb effect in the collision 



of lighter nuclei the effect of incident energy is almost same as one move from low energy to 



8 



high energy. While in case of heaviest nuclei (Au+Au) relative ratio of fragment production 
at low energy is FN: LMF: IMF :: 9: 3: 1 and at high energy FN: LMF: IMF :: 83: 9: 
1. There is large change in the formation of fragment, due to large compression at high 
energy and Coulomb repulsion. We study that the rate of change of fragment production 



is decrease as we go beyond 400 MeV/nucleon. It is also experimentally o 



Dserved that a 



little change in fragmentation yield takes place beyond 400 MeV/nucleon [l7j. The slope of 
fragment production in all the cases below below saturation time is high due to the Pauli 
blocking. 

On considering the influence of symmetry energy on fragment production, we analyse FN 
for neutron rich nuclei is almost independent from the affect of symmetry energy. But 
one can see a measureable effect of symmetry energy for neutron deficient nuclei at 100 
MeV/nucleon. In LMF production influence of symmetry energy is large as compared to 
FN. We got pronounced effect of different symmetry energies in both the cases as shown in 
second horizontal panel. On observing IMF's production for light nuclei we can say IMF 
multiplicity increases with increase in the stiffness of symmetry energy. But for heavy nuclei 
production of IMF is independent of symmetry energy and we observed measureable effect of 
symmetry energy only at low energy. While for light nuclei we see the influence of symmetry 
energy at all energy range as shown in figure. 

So far we have seen the effect of incident energy, symmetry energy on fragment production. 
To see the effect of symmetry in more effective way we include the isoscaling parameter 
and cross- sect ion. Here we fix cross-section at 20 and 55 mb to study the effect of multi- 
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FIG. 3: Effect of cross-section with incident energy on production of free nucleon at different 
symmetry energy between a neutron rich nuclei Au and neutron deficient nuclei Ne. 

fragmentation. For comparison, we shall also use an isospin dependent cross-section where 
cr„p is more as compare to app, Onn- We see the effect of cross-section on production of FN 
and LMF for two different multifragmentation reactions i.e. ratio of production of FN and 
LMF studied between a neutron rich nuclei Au and a neutron deficient nuclei Ne for incident 
energies 50, 100, 200, 400, 600, 1000 MeV/nucleon for different symmetry energy. 
In first panel we consider isospin dependent cross- sect ion. In second and third panel we use 
constant isospin independent cross-section with value 20 mb and 55 mb respectively. One 
can see that fragment production in addition to symmetry energy is affected by different 
cross-section. 

In the figure 3 and 4 we display the relative production of FN is more than LMF's in all 
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FIG. 4: Effect of cross-section with incident energy on production of light mass fragments at 
different symmetry energy between a neutron rich nuclei Au and neutron deficient nuclei Ne. 

the cases of isospin dependent and fixed cross-section. We can justify our study on the basis 



of the Cugnon parameterization 221]. For the energy range 50 to 1000 MeV/nucleon, elastic 
collisions take place. In figure 2 we see the influence of symmetry energy on production of 
FN and LMF, but when we include the cross section, the slope of relative ratio of neutron 
rich and deficient nuclei decrease with increase in incident energy as shown in fig. 3 and 
4. Moreover the role of cross-section decreases with increase in energy. The behaviour of 
relative yield curve at 7 = 2 shift due to the energy effect, the relative yield ratio at aiso 
and (755 is approximately 10 at high energy, while at a2o its value decrese for FN. For LMF 
the relative ratio at aiso and (755 is approximately 5 at low energy and 8 for high energy. 
In figure 5 we compare our study of IMF's production at different energies of Au + Au 
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collision with experimental data of ALADIN for central collisions |l7l . [iSj. The production 
of IMF's decrease with incident energy, because at high energy the nuclear matter shatter 
down in small pieces which leads to more production of FN and LMF at the cost of decrease 
in the production of IMF's. In figure 5 we have displayed the results at cross-section 20 
mb, 55 mb, and aiso in the absence (7 = 0) and presence (7 = 0.66) of density dependence 
of symmetry energy. One can see that although the theoretical results are in not in close 
agreement with the data, but the trend of result in similar to experimental finding. The 
mismatch in the trend may be due to the binding energy and momentum parameter. As 
stated earlier cluster and fragments in this study were generated by simple minimum span- 
ning tree method. This method binds two nucleons based on spatial correlation regardless 
of the fact whether proper binding energy is achieved or not. The new methods based on 



binding energy creteria can shed light on this aspect 23|. It was noted by many authors 



24j | that the various forms of equation of state as well as momentum dependent interactions 



can also make different impact on the results. 



IV. SUMMARY 

In this paper we have used isospin dependent quantum molecular dynamics model 
to study the influence of symmetry energy and cross-section on multifragmentation at 
neutron rich (Au) nuclei as well as neutron deficient (Ne) nuclei at incident energy 50-1000 
MeV/nucleon. All the simulations carried at central mass frame with scaled impact 
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FIG. 5: Comparison of average multiplicity of intermediate mass fragments (IMFs) with ALADIN 
data at incident energies of 100, 200, 400, 600, 1000 MeV/nucleon as a function of energy. 

parameter 0.3. We used momentum independent and soft equation of state. Our study 
showed that the multiphcity of fragments depends upon incident energy, symmetry energy 
influences the fragment production, effect of symmetry energy is more for neutron deficient 
nuclei than the neutron rich nuclei and cross section affects the fragment production. 
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